DESY 94-044 
March 1994 



0(a 2 L 2 ) radiative corrections to deep inelastic 
ep scattering for different kinematical 

variables 



Johannes Bliimlein 



DESY-Institut fur Hochenergiephysik, Zeuthen, 
Platanenallee 6, D-15735 Zeuthen, Germany 



Abstract 

The QED radiative corrections are calculated in the leading log approximation up to 
0{a 2 ) for different definitions of the kinematical variables using jet measurement, the 
'mixed' variables, the double angle method, and a measurement based on 6 e and yjB- 
Higher order contributions due to exponentiation of soft radiation are included. 



1 Introduction 



The measurement of deep inelastic ep scattering at HERA will allow to extend the kinematical 
range to x ~ 1(T 4 and Q 2 ~ C(10 4 ) GeV 2 §. Because of differences in the detector response 
and resolution 0] different ways to measure the kinematical variables x, y and Q 2 have to be 
used to cover the full kinematical range. Furthermore, owing to the different properties of de- 
tectors the experiments HI and ZEUS || apply different methods in measuring the kinematical 
variables. 

The QED radiative corrections in 0(a) differ significantly for the various choices of kine- 
matical variables as shown in || . This is due to the full (or partial) integration of the radiated 
photon's phase space for different kinematical situations according to the experimental require- 
ments. As shown in explicit comparisons between complete and leading log calculations in 0(a) 
the results agree to the per cent level for the case of lepton measurement |5], |(J, jet measurement 
0, H, and the case of mixed variables |7|, || 0. In part of the phase space the 0(a) corrections 
are sizeable. Therefore, it is necessary to calculate the 2nd order corrections to check whether 
the results obtained are numerically stable. 

In this letter the 0(a 2 ) contributions are calculated in the leading log approximation for the 
case of jet measurement, mixed variables, the double angle method JTTL and a measurement 
based on 6 e and yjacquet-Biondei = Vjb & The latter two methods have been applied by the ZEUS 
collaboration recently in measuring the deep inelastic neutral current cross section JTB|. Since 
the 0(aL) radiative corrections for these cases have not been given before the corresponding 
numerical results will be included also. Besides of the second order leading log corrections 
higher order corrections due to soft exponentiation are given. 



2 0(a 2 L 2 ) corrections 

The relevant contributions arise from the set of all possible collinear configurations in the 
diagrams up to 0(a 2 ). They are associated with the 'massless' fermion lines (m 2 ,m 2 -C Q 2 ). 
One may group the different processes into: 

1. bremsstrahlung (diagrams la,b, figure 1) 

2. e + e~ pair creation (diagram lc), and 

3. // pair production, / = l~ ,u,d, s, c, b, (diagram Id). 

Among the possible configurations contributions due to initial state and final state radiation 
both from the lepton and quark lines emerge. Furthermore, the Compton scattering of a nearly 
real intermediary photon, emitted from a quark or proton line and scattered off the electron, 
|| may be considered as a part of the radiative corrections. Since the experimental signature 
of the latter process consists of two nearly balanced electromagnetic showers accompanied by 
low hadronic activity near the beam pipe, these events are generally not included into the 
deep inelastic sample. On the other hand, they may be used to measure the proton structure 

1 Also the case of 'hadronic' variables was investigated ^, [T(|. However, these variables are not accessible in 
deep inelastic scattering experiments. 

2 For leptonic variables these corrections have been calculated in JL2J. 
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functions at small x and Q 2 JT4J. Because of the finite detector resolutions collinear final state 
radiation is hardly resolved and may be treated integrally [15 1 in the experimental measurement. 
Due to this also terms with both collinear initial and final state radiation do not emerge. The 



initial state radiation from the quark line can be summed using Altarelli-Parisi equations [16 



with kernels modified due to the 1st and 2nd order QED corrections, which yields a modification 
of the parton distributions at the level of 1% i. 

Because of this, we will consider only the effect of initial state radiation in the following. 
For the scattering cross section up to 0(a 2 ) one thus obtains: 
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using standard techniques known from QCD □. Here, d 2 a^ /dxdy denotes the Born cross section 
for deep inelastic neutral or charged current ep reactions! (cf. e.g. (3,4) in |7|]), and the Jacobian 
J is given by 

dx/dx dy/dx 
dx/dy dy/dy 



J(x,y,z) 



(2) 



The shifted variables x = Q 2 /ys, y, Q 2 , s, and the lower bounds of z, Zq, are listed in table 1 
for the different choices of measurement. 

Both for the double angle method and the measurement based on 9 e and yjB, the rescaled 
variables x and Q 2 vanish for z — > Zq (cf. table 1). If no futher cut is imposed this means that 
for both methods the radiative corrections are strongly influenced by the non-perturbative 
range x, Q 2 — > in d 2 of dxdy. In studies by the ZEUS collaboration, where these two methods 
have been applied, however, the cut 



2E P 



E'Jl - cos6 e ) + Ej(l - cosOj) >A 



(3) 



on the reconstructed electron beam energy was imposed, choosing ^4 = 35 GeV. Here, # e j 
denote the scattering angle of the electron and the jet angle, respectively, and E e , E' e , Ej are 
the energies of the initial electron, scattered electron, and the jet in the laboratory frame. This 
cut yields z = A/{2E e ) in the case of the double angle method and z = ma.x{A/ (2E e ), y} 



3 Note, that the current uncertainty in the parton distributions is larger than this in most of the kinematical 
range 

4 These methods have also been used for the calculation of QED corrections to e + e annihilation processes 
successfully pl| . 

5 Note, that QED loop effects, as e.g. the polarization of the intermediary boson, are not included in ([I]). 
One may easily account for them using 'dressed' Born cross sections instead of da^ /dxdy. 
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for the measurement based on 6 e and yjB- For the double angle method the cut excludes the 
contributions of the range x, Q 2 — > globally, whereas they are still present in the case of the 
0e — Vjb measurement for A/(2E e ) ~ y. 
The splitting functions in (H) are: 
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Here, P« and P« are the splitting functions of a photon into an electron, and vice versa. 
® denotes the Mellin convolution, ef is the fermion charge, and N c (f) = 3 for quarks, 
N c (f) = 1 for leptons, respectively. (|7|) is obtained from an expansion of otQED(s) calcu- 
lating the polarization operator in the on-mass-shell scheme for leptons and quarks, which 
can not be done reliably in perturbative QCD in the case of the light quarks. However, 
one may try to parametrize C(qed{s, wi/J in terms of effective quark masses from a fit to 
R(s) = cr(e + e~ — > hadrons)/<r(e + e~ — > // + fi~) as done in f22[ 
tion here and use the parameters m u = 62 MeV, = 83 MeV, m s 
and rrib = 4.5 GeV as obtained in |22| in the subsequent analysis. 



We will apply this descrip- 
215 MeV, m r = 1.5 GeV, 



The soft exponentiation beyond the second order can be calculated using the solution of 
the non-singlet evolution equation for the electron or positron distribution in the soft range 
(z — > 1) for running Oqed- It is given by |23| 



D NS (z,Q 2 ) =C(1 
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(9) 



Expanding (||) up to 0(a 2 ), one obtains terms in the limit z — > 1 which are contained in (|l|) 
already. Apart from bremsstrahlung terms also the e + e~ pair production term in ([!]) is included 
for m e — > i. To account for the terms of 0(a 3 ) and higher orders contained in we use 



f?, 2 ' *°"feC 2 ) = D KS (z,Q' 
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6 One might generalize (^) using also the heavier fermions u,d, fx, s, ... in the running coupling constant. 
However, apart from of the mass threshold in (J7|) the smaller logarithm associated with these terms will yield 
an even smaller contribution. 
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anJ 

= / dzP£ 2 \z,Q 2 ) \0 {Z - Zo )J( X ,y, z) — - _) (11) 

Beginning with the second order corrections in a terms due to e~ — > e + conversion in the 
initial state contribute. They result from diagram lc in 0(a 2 ) where now the outgoing fermion 
collinear with the incoming electron is a positron. The conversion rate equals 

P(z,QV-^e+)=(^) 2 l„'(g)pM W (12) 

This quantity is illustrated in figure 2 numerically. For Q 2 ~ lOGeV 2 , P(z,Q 2 ;e~ — > e + ) 
may reach values of 0(1) for z ~ 10~ 4 . Depending on the detector resolution the charge 
of the scattered lepton in deep inelastic ep collisions might only be measurable in part of 
the kinematical range. If the charge of the recorded final state electron is ignored in the 
measurement, also the contribution 

d 2 a (2,e--,e+) fl ^2^(0) 

dzP(z, Q 2 - e -> e + )J(x, y, z) \ x=s . $ s (13) 



dxdy Jz " dxdy 

has to be included in the radiative corrections. The corresponding contributions are given for 
the different cases of measurement separately below. 



3 Numerical results 



In the numerical calculation the quark distributions are parametrized in the DIS scheme using 
the MRS D~ parton distributions [17] to illustrate the correction functions 5i(x,y) for the 
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different type of measurements. Equivalent results are obtained using the recent MRSH 
CTEQ H, or GRV |0| distributions. 

The value of Q 2 for the structure functions in the radiative orders in (|1|) may become smaller 
than the usual lower bound Ql ~ 4 GeV 2 for which the chosen parton parametrization applies. 
Because of Lorentz invariance the structure functions F% 3 (x, Q 2 ) vanish with Q 2 — > 0. Various 



parametrizations of this behaviour have been proposed recently f24 |. Here, for an illustration, 
we will adopt a multiplicative factor 



x 



1 -exp(-A 2 Q 2 



with A = 3.37 GeV 



-2 



(14) 



for all structure functions in the numerical calculations. 

Figure 3 illustrates the 0(a 2 L 2 ) corrections for the case of jet measurement for neutral 
current deep inelastic e~p scattering. The corrections are positive and grow with x and y. 
Compared with the 1st order results they diminish the negative 0{a) results in size. In the 
low x range x ~ 0.01 the 2nd order corrections are found to be less than 1 %. Still smaller 
values are found for the 2nd order contribution due to e~ — > e + conversion, which amount to 
relative corrections < 4 • 10~ 4 only. 



7 To obtain consistency between eqs. ( jio| ) and the limit m e — > in (0) is taken. 
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The 2nd order corrections for jet measurement in the case of charged current deep inelastic 
e~p scattering are the same as those of the corresponding neutral current reactions for large 
x, because this range is dominated by the soft and virtual contributions, cf. figure 4. For 
smaller values of x ~ 0.01 they are smaller than 1 % and may become negative, as also the 1st 
order corrections change sign for small x J/J. The contributions due to e~ — > e + conversion are 
somewhat larger than in the case of neutral current reactions. 

For the measurement based on mixed variables, i.e. Q 2 measured at the leptonic and yjs at 
the hadronic vertex (see figure 5), the 2nd order corrections to the neutral current cross section 
turn out to be less than 1 % in the range y ~ 0.2. They become smaller with descending x 
values. For x ~ 0.5 they reach 10 % in the range of small y. Note that contrary to the jet 
measurement the contributions due to e~ — > e + conversion become of comparable order as the 
2nd order e~ — *> e~ corrections for small y in the case of mixed variables. 

The 0{a) and 0(a 2 ) corrections for the neutral current deep inelastic scattering cross 
section for a measurement using the double angle method are shown in figure 6. Both the 1st 
and 2nd order corrections are nearly flat in y. The corrections become smaller with descending 
values of x and are smaller than 3 % for x ~ 0.01. Due to these properties the double angle 
method appears to be an ideal choice to measure the neutral current cross section from the 
point of view of radiative corrections. The contributions from e~ — > e + conversion are less than 
2 ■ 10~ 4 . 

Figure 7 shows the O(o) and 0(a 2 ) corrections for the neutral current deep inelastic scat- 
tering cross section, where the scattering angle 9 e of the electron ond yjB are used to define the 
kinematical variables. In this case the corrections grow rapidly in the range y m Aj (2E e ) and 
become larger than 100 % for y ~ A/(2E e ) already for x = 0.01. Even larger corrections are 
obtained for higher x values. This is caused by the shift of the rescaled variables x — > and 
Q 2 — > in the vicinity of z m y. The relative correction becomes smaller for very low values of 
x, since the differential cross sections at the rescaled variables and the non-rescaled variables 
become less different in value there. The 2nd order corrections are negative. The contributions 
due to — ► e + conversion grow with rising values of x and amount to 0(2%) of the Born cross 
section for x ~ 0.5 and larger values of y. Because in part of the kinematical range the correc- 
tion 8nc{x, y) is significantly influenced by the behaviour of da^ jdxdy in the non-perturbative 
range, a reliable calculation of 5rfc(%,y) is only possible knowing the behaviour da^ jdxdy in 
the range of Q 2 ,x — ► 0. Due to this the plot of 5nc{ x ,v) m figure 7 serves only to show the 
principal behaviour under some assumptions made for da^ jdxdy. If, on the other hand, the 
cross section da/dxdy, which includes the radiative corrections, is well understood with respect 
to the experimental systematics, one may later try to use the measurement based on 8 e and 
yjB as a method to unfold da^ jdxdy in the non-perturbative range. 

In summary, we found that the 0(a 2 L 2 ) corrections, supplemented by the soft exponentia- 
tion, for the neutral and charged current deep inelastic ep scattering cross sections amount to 
~ 10% for the cases of jet measurement, mixed variables, and the double angle method. These 
contributions are of the order of the difference between the complete 0(a) and the leading log 
result 0(aL). In the case of the measurement based on 9 e and yjB, contributions from the 
non-perturbative range x, Q 2 — > influence both the 1st and 2nd order corrections. 
Acknowledgement. I would like to thank Paul Soding for reading the manuscript. 
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Table 1: Description of parameters in eq. ([!]) for different types of cross section measurement 
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Figure 1: Diagrams contributing to the radiative corrections up to 0(a 2 L 2 ^ 
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Figure 3: Leptonic initial state radiative corrections Snc( x jU) = (da^ 2+>2 ' so ^ / dxdy) / (da / dxdy) in 
LLA for e~p deep inelastic scattering in the case of jet measurement for y/s = 314 GeV, .A = 0, and 
Q 2 > 5 GeV 2 . Full lines: 0(a 2 ) corrections. Dashed lines: contributions due to e~ — ► e + conversion 
eq. (PI), 5 e N( ^* e+ (x, y) = (da^ 2,e ~* e ) /dxdy) /{da /dxdy) scaled by x50; upper line: x = 0.01, middle 
line: x = 0.0001, lower line x = 0.9. 
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Figure 4: 5cc(x,y) = (da^ >2 ' so ^^ / dxdy) / {da^u / dxdy) for deep inelastic e~p scattering in the case 
of jet measurement. Dotted lines: ^c^ 6 ( x iV)- The other parameters are the same as in figure 3. 
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Figure 5: 5nc{x,v) for the case of mixed variables. Dotted lines: &%c* e (x,y); upper line: x = 0.5, 
lower line x = 0.01. The other parameters are the same as in figure 3. 
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Figure 6: 5nc(x,v) for the case of the double angle method for A = 35 GeV. 
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as in figure 3. 



0.01, lower line: x = 0.0001. The other parameters are the same 
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Figure 7: 5nc(x,v) for the measurement based on 9 e and yjs for A = 35 GeV. Full lines: 

Dashed lines: 8%^ e+ (x,y); upper line: x = 0.5, 
The other parameters are the same as in figure 3. 
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